A novel homozygous splice site variant in ARL2BP causes a syndromic autosomal recessive rod-cone dystrophy with situs inversus, asthenozoospermia, unilateral renal agenesis and microcysts.

BACKGROUND: This report presents a clinical case of syndromic rod-cone dystrophy due to a splice site variant in the ARL2BP gene causing situs inversus, asthenozoospermia, unilateral renal agenesis and microcysts. The presence of renal agenesis and cryptorchidism expands the clinical manifestations due to ARL2BP variants. The detailed, long-term follow-up contributes valuable insights into disease progression, aiding clinical diagnosis and patient management. CASE PRESENTATION: The male patient complained of photophobia as the first symptom when he was 20 years old followed by nyctalopia, loss of central visual acuity and peripheral visual field ten years later. Genetic analysis identified a likely pathogenic homozygous variant (c.294-1G > C) involving the splicing acceptor site of intron 4. Reported symptoms together with full-field stimulus threshold testing, electroretinogram and advanced multimodal imaging allowed us to recognize the typical characteristics of a mixed retinal dystrophy. Despite the end-stage retinal disease, this patient still retained a useful residual vision at 63 years and had a slow disease progression during the last 5 years of evaluation. DISCUSSION AND CONCLUSIONS: Our findings underscore the variable clinical presentation of ARL2BP variants, emphasizing the importance of a nuanced approach in diagnosing and managing patients. The presence of renal cysts warrants consideration of a differential diagnosis, particularly with Senior-Loken (SLS), Bardet-Biedl (BBS) and Joubert syndromes (JS) but also with Short Rib Thoracic Dysplasia 9, highlighting the need for careful phenotypic evaluation in these cases.

Spliceosomal helicases DDX41/SACY-1 and PRP22/MOG-5 both contribute to proofreading against proximal 3' splice site usage.

RNA helicases drive necessary rearrangements and ensure fidelity during the pre-mRNA splicing cycle. DEAD-box helicase DDX41 has been linked to human disease and has recently been shown to interact with DEAH-box helicase PRP22 in the spliceosomal C* complex, yet its function in splicing remains unknown. Depletion of DDX41 homolog SACY-1 from somatic cells has been previously shown to lead to changes in alternative 3' splice site (3'ss) usage. Here, we show by transcriptomic analysis of published and novel data sets that SACY-1 perturbation causes a previously unreported pattern in alternative 3' splicing in introns with pairs of 3' splice sites ≤18 nt away from each other. We find that both SACY-1 depletion and the allele sacy-1(G533R) lead to a striking unidirectional increase in the usage of the proximal (upstream) 3'ss. We previously discovered a similar alternative splicing pattern between germline tissue and somatic tissue, in which there is a unidirectional increase in proximal 3'ss usage in the germline for ∼200 events; many of the somatic SACY-1 alternative 3' splicing events overlap with these developmentally regulated events. We generated targeted mutant alleles of the Caenorhabditis elegans homolog of PRP22, mog-5, in the region of MOG-5 that is predicted to interact with SACY-1 based on the human C* structure. These viable alleles, and a mimic of the myelodysplastic syndrome-associated allele DDX41(R525H), all promote the usage of proximal alternative adjacent 3' splice sites. We show that PRP22/MOG-5 and DDX41/SACY-1 have overlapping roles in proofreading the 3'ss.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs.

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames and without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs may account for proteins that have been attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

Estimating the proportion of nonsense variants undergoing the newly described phenomenon of manufactured splice rescue.

A recent report described a nonsense variant simultaneously creating a donor splice site, resulting in a truncated but functional protein. To explore the generalizability of this unique mechanism, we annotated >115,000 nonsense variants using SpliceAI. Between 0.61% (donor gain delta score >0.8, for high precision) and 2.57% (>0.2, for high sensitivity) of nonsense variants were predicted to create new donor splice sites at or upstream of the stop codon. These variants were less likely than other nonsense variants in the same genes to be classified as pathogenic/likely pathogenic in ClinVar (p < 0.001). Up to 1 in 175 nonsense variants were predicted to result in small in-frame deletions and loss-of-function evasion through this "manufactured splice rescue" mechanism. We urge caution when interpreting nonsense variants where manufactured splice rescue is a strong possibility and correlation with phenotype is challenging, as will often be the case with secondary findings and newborn genomic screening programs.

Systematic Minigene-Based Splicing Analysis and Tentative Clinical Classification of 52 CHEK2 Splice-Site Variants.

BACKGROUND: Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes. METHODS: A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells. RESULTS: Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition. CONCLUSIONS: Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants.

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Impact on splicing in Saccharomyces cerevisiae of random 50-base sequences inserted into an intron.

Intron splicing is a key regulatory step in gene expression in eukaryotes. Three sequence elements required for splicing-5' and 3' splice sites and a branchpoint-are especially well-characterized in Saccharomyces cerevisiae, but our understanding of additional intron features that impact splicing in this organism is incomplete, due largely to its small number of introns. To overcome this limitation, we constructed a library in S. cerevisiae of random 50-nt (N50) elements individually inserted into the intron of a reporter gene and quantified canonical splicing and the use of cryptic splice sites by sequencing analysis. More than 70% of approximately 140,000 N50 elements reduced splicing by at least 20%. N50 features, including higher GC content, presence of GU repeats, and stronger predicted secondary structure of its pre-mRNA, correlated with reduced splicing efficiency. A likely basis for the reduced splicing of such a large proportion of variants is the formation of RNA structures that pair N50 bases-such as the GU repeats-with other bases specifically within the reporter pre-mRNA analyzed. However, multiple models were unable to explain more than a small fraction of the variance in splicing efficiency across the library, suggesting that complex nonlinear interactions in RNA structures are not accurately captured by RNA structure prediction methods. Our results imply that the specific context of a pre-mRNA may determine the bases allowable in an intron to prevent secondary structures that reduce splicing. This large data set can serve as a resource for further exploration of splicing mechanisms.

Novel CWF19L1 mutations in patients with spinocerebellar ataxia, autosomal recessive 17.

Spinocerebellar ataxia, autosomal recessive-17 (SCAR17) is a rare hereditary ataxia characterized by ataxic gait, cerebellar signs and occasionally accompanied by intellectual disability and seizures. Pathogenic mutations in the CWF19L1 gene that code for CWF19 like cell cycle control factor 1 cause SCAR17. We report here two unrelated families with the clinical characteristics of global developmental delay, cerebellar ataxia, pyramidal signs, and seizures. Cerebellar atrophy, and T2/FLAIR hypointense transverse pontine stripes were observed in brain imaging. Exome sequencing identified novel homozygous mutations including a splice acceptor site variant c.1375-2 A > G on intron 12 in a male patient and a single nucleotide variant c.452 T > G on exon 5 resulting in a missense variant p.Ile151Ser in the female patient from two unrelated families, respectively. Sanger sequencing confirmed the segregation of these variants in the family members with autosomal recessive inheritance. Transcript analysis of the splice site variant revealed activation of a novel cryptic splice acceptor site on exon 13 resulting in an alternative transcription with a loss of nine nucleotides on exon 13. Translation of this transcript predicted an in-frame deletion of three amino acids p.(459_461del). We also observed a novel exon 13 skipping which results in premature termination of the protein product. Our study expands the phenotype, radiological features, and genotypes known in SCAR17.

Sequence variations affect the 5' splice site selection of plant introns.

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Pathogenicity analysis and splicing rescue of a classical splice site variant (c.1343+1G>T) of CNOT1 gene associated with neurodevelopmental disorders.

Mutations in the CNOT1 gene lead to an incurable rare neurological disorder mainly manifested as a clinical spectrum of intellectual disability, developmental delay, seizures, and behavioral problems. In this study, we investigated a classical splice site variant of CNOT1 (c.1343+1G>T) associated with neurodevelopmental disorders, which was a master regulator, orchestrating gene expression, RNA deadenylation, and protein ubiquitination. To link CNOT1 dysfunction with the neurodevelopmental phenotype observed in a patient, in vitro minigene assay was used to verify the effect of CNOT1 gene splice site variant c.1343+1G>T on mRNA splicing. We also explored the impact of transient transfection introducing modified U1 snRNA on correcting the splicing variant. Through minigene expression in mammalian cells, we demonstrated that the variant induced complete exon 12 skipping, which explained the patient's clinical condition and provided additional genetic diagnosis evidence for the clinical significance of the variant. Moreover, we confirmed that the aberrant splice pattern could be partially corrected by the modified U1 snRNA at the mRNA level, which provided strong evidence for the therapeutic potential of modified U1 snRNA in neutralizing the hazardous effect of incorrect splicing patterns.

Mouse nuclear RNAi-defective 2 promotes splicing of weak 5' splice sites.

Removal of introns during pre-mRNA splicing, which is central to gene expression, initiates by base pairing of U1 snRNA with a 5' splice site (5'SS). In mammals, many introns contain weak 5'SSs that are not efficiently recognized by the canonical U1 snRNP, suggesting alternative mechanisms exist. Here, we develop a cross-linking immunoprecipitation coupled to a high-throughput sequencing method, BCLIP-seq, to identify NRDE2 (nuclear RNAi-defective 2), and CCDC174 (coiled-coil domain-containing 174) as novel RNA-binding proteins in mouse ES cells that associate with U1 snRNA and 5'SSs. Both proteins bind directly to U1 snRNA independently of canonical U1 snRNP-specific proteins, and they are required for the selection and effective processing of weak 5'SSs. Our results reveal that mammalian cells use noncanonical splicing factors bound directly to U1 snRNA to effectively select suboptimal 5'SS sequences in hundreds of genes, promoting proper splice site choice, and accurate pre-mRNA splicing.

Recursive splicing discovery using lariats in total RNA sequencing.

Recursive splicing is a non-canonical splicing mechanism in which an intron is removed in segments via multiple splicing reactions. Relatively few recursive splice sites have been identified with high confidence in human introns, and more comprehensive analyses are needed to better characterize where recursive splicing happens and whether or not it has a regulatory function. In this study, we use an unbiased approach using intron lariats to search for recursive splice sites in constitutive introns and alternative exons in the human transcriptome. We find evidence for recursive splicing in a broader range of intron sizes than previously reported and detail a new location for recursive splicing at the distal ends of cassette exons. In addition, we identify evidence for the conservation of these recursive splice sites among higher vertebrates and the use of these sites to influence alternative exon exclusion. Together, our data demonstrate the prevalence of recursive splicing and its potential influence on gene expression through alternatively spliced isoforms.

Generation of tandem alternative splice acceptor sites and CLTC haploinsufficiency: A cause of CLTC-related disorder.

Tandem splice acceptors (NAGNn AG) are a common mechanism of alternative splicing, but variants that are likely to generate or to disrupt tandem splice sites have rarely been reported as disease causing. We identify a pathogenic intron 23 CLTC variant (NM_004859.4:c.[3766-13_3766-5del];[=]) in a propositus with intellectual disability and behavioral problems. By RNAseq analysis of peripheral blood mRNA, this variant generates transcripts using cryptic proximal splice acceptors (NM_004859.4: r.3765_3766insTTCACAGAAAGGAACTAG, and NM_004859.4:r.3765_3766insAAAGGAACTAG). Given that the propositus expresses 38% the level of CLTC transcripts as unaffected controls, these variant transcripts, which encode premature termination codons, likely undergo nonsense mediated mRNA decay (NMD). This is the first functional evidence for CLTC haploinsufficiency as a cause of CLTC-related disorder and the first evidence that the generation of tandem alternative splice sites causes CLTC-related disorder. We suggest that variants creating tandem alternative splice sites are an underreported disease mechanism and that transcriptome-level analysis should be routinely pursued to define the pathogenicity of such variants.

All reported non-canonical splice site variants in GLA cause aberrant splicing.

BACKGROUND: Fabry disease is an X-linked lysosomal storage disorder caused by insufficient α-galactosidase A (GLA) activity resulting from variants in the GLA gene, which leads to glycosphingolipid accumulation and life-threatening, multi-organ complications. Approximately 50 variants have been reported that cause splicing abnormalities in GLA. Most were found within canonical splice sites, which are highly conserved GT and AG splice acceptor and donor dinucleotides, whereas one-third were located outside canonical splice sites, making it difficult to interpret their pathogenicity. In this study, we aimed to investigate the genetic pathogenicity of variants located in non-canonical splice sites within the GLA gene. METHODS: 13 variants, including four deep intronic variants, were selected from the Human Gene Variant Database Professional. We performed an in vitro splicing assay to identify splicing abnormalities in the variants. RESULTS: All candidate non-canonical splice site variants in GLA caused aberrant splicing. Additionally, all but one variant was protein-truncating. The four deep intronic variants generated abnormal transcripts, including a cryptic exon, as well as normal transcripts, with the proportion of each differing in a cell-specific manner. CONCLUSIONS: Validation of splicing effects using an in vitro splicing assay is useful for confirming pathogenicity and determining associations with clinical phenotypes.

A deep intronic variant in DNM1 in a patient with developmental and epileptic encephalopathy creates a splice acceptor site and affects only transcript variants including exon 10a.

DNM1 developmental and epileptic encephalopathy (DEE) is characterized by severe to profound intellectual disability, hypotonia, movement disorder, and refractory epilepsy, typically presenting with infantile spasms. Most of the affected individuals had de novo missense variants in DNM1. DNM1 undergoes alternative splicing that results in expression of six different transcript variants. One alternatively spliced region affects the tandemly arranged exons 10a and 10b, producing isoforms DNM1A and DNM1B, respectively. Pathogenic variants in the DNM1 coding region affect all transcript variants. Recently, a de novo DNM1 NM_001288739.1:c.1197-8G > A variant located in intron 9 has been reported in several unrelated individuals with DEE that causes in-frame insertion of two amino acids and leads to disease through a dominant-negative mechanism. We report on a patient with DEE and a de novo DNM1 variant NM_001288739.2:c.1197-46C > G in intron 9, upstream of exon 10a. By RT-PCR and Sanger sequencing using fibroblast-derived cDNA of the patient, we identified aberrantly spliced DNM1 mRNAs with exon 9 spliced to the last 45 nucleotides of intron 9 followed by exon 10a (NM_001288739.2:r.1196_1197ins[1197-1_1197-45]). The encoded DNM1A mutant is predicted to contain 15 novel amino acids between Ile398 and Arg399 [NP_001275668.1:p.(Ile398_Arg399ins15)] and likely functions in a dominant-negative manner, similar to other DNM1 mutants. Our data confirm the importance of the DNM1 isoform A for normal human brain function that is underscored by previously reported predominant expression of DMN1A transcripts in pediatric brain, functional differences of the mouse Dnm1a and Dnm1b isoforms, and the Dnm1 fitful mouse, an epilepsy mouse model.

A deep intronic TCTN2 variant activating a cryptic exon predicted by SpliceRover in a patient with Joubert syndrome.

The recent introduction of genome sequencing in genetic analysis has led to the identification of pathogenic variants located in deep introns. Recently, several new tools have emerged to predict the impact of variants on splicing. Here, we present a Japanese boy of Joubert syndrome with biallelic TCTN2 variants. Exome sequencing identified only a heterozygous maternal nonsense TCTN2 variant (NM_024809.5:c.916C >T, p.(Gln306Ter)). Subsequent genome sequencing identified a deep intronic variant (c.1033+423G>A) inherited from his father. The machine learning algorithms SpliceAI, Squirls, and Pangolin were unable to predict alterations in splicing by the c.1033+423G>A variant. SpliceRover, a tool for splice site prediction using FASTA sequence, was able to detect a cryptic exon which was 85-bp away from the variant and within the inverted Alu sequence while SpliceRover scores for these splice sites showed slight increase (donor) or decrease (acceptor) between the reference and mutant sequences. RNA sequencing and RT-PCR using urinary cells confirmed inclusion of the cryptic exon. The patient showed major symptoms of TCTN2-related disorders such as developmental delay, dysmorphic facial features and polydactyly. He also showed uncommon features such as retinal dystrophy, exotropia, abnormal pattern of respiration, and periventricular heterotopia, confirming these as one of features of TCTN2-related disorders. Our study highlights usefulness of genome sequencing and RNA sequencing using urinary cells for molecular diagnosis of genetic disorders and suggests that database of cryptic splice sites predicted in introns by SpliceRover using the reference sequences can be helpful in extracting candidate variants from large numbers of intronic variants in genome sequencing.

Systematic analysis of CNGA3 splice variants identifies different mechanisms of aberrant splicing.

Achromatopsia is an autosomal recessive cone photoreceptor disease that is frequently caused by pathogenic variants in the CNGA3 gene. Here, we present a systematic functional analysis of 20 CNGA3 splice site variants detected in our large cohort of achromatopsia patients and/or listed in common variant databases. All variants were analyzed by functional splice assays based on the pSPL3 exon trapping vector. We demonstrated that ten variants, both at canonical and non-canonical splice sites, induced aberrant splicing, including intronic nucleotide retention, exonic nucleotide deletion and exon skipping, resulting in 21 different aberrant transcripts. Of these, eleven were predicted to introduce a premature termination codon. The pathogenicity of all variants was assessed based on established guidelines for variant classification. Incorporation of the results of our functional analyses enabled re-classification of 75% of variants previously classified as variants of uncertain significance into either likely benign or likely pathogenic. Our study is the first in which a systematic characterization of putative CNGA3 splice variants has been performed. We demonstrated the utility of pSPL3 based minigene assays in the effective assessment of putative splice variants. Our findings improve the diagnosis of achromatopsia patients, who may thus benefit from future gene-based therapeutic strategies.